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Введение. Цилиндрические и сферические оболочки широко используются в технике. Они подвергаются 
внутреннему и/или внешнему давлению и тепловому воздействию. Распределение напряжений и деформаций в 
упругопластических оболочках изучалось многими исследователями. Большинство работ связано с 
использованием условий пластичности Мизеса, максимального касательного, максимального приведённого 
напряжения. Эти условия не учитывают зависимость от первого инварианта тензора напряжений и знака 
третьего инварианта девиатора напряжений. Для тел со сферической и цилиндрической симметрией при 
осесимметричном тепловом и силовом воздействии в ряде случаев удается получить численно-аналитические 
решения для напряжений, перемещений и деформаций. 

Материалы и методы. Решение задачи о состоянии толстостенной упругопластической оболочки проводится в 
рамках теории малых деформаций. Предложено условие пластичности, учитывающее зависимость от трех 
инвариантов тензора напряжений, а также знак третьего инварианта девиатора напряжений и трансляционное 
упрочнение материала. Решается несвязная термоупругопластическая задача. Для оценки напряжений в 
области упругого состояния сферической оболочки вводится эквивалентное напряжение, равное выбираемой 
функции пластичности. В качестве метода верификации напряженного состояния используется построение 
годографа вектора напряжений. 

Результаты исследования. Для линейных функций пластичности задача имеет аналитическое решение. 
Получено решение, учитывающее упрочнение материала. Определены аналитические и графические 
зависимости между параметрами внешнего воздействия для упругого и упругопластического состояния шара. В 
случае комбинированной нагрузки возможны варианты, когда пластическая область зарождается на 
внутренней, внешней границах шара или между этими границами. 

Обсуждение и заключения. Результаты вычислений показали, что учет пластической сжимаемости и 
зависимости предела пластичности от температуры может оказать существенное влияние на напряженное и 
деформированное состояние полого шара. При этом учет первого инварианта тензора напряжений в условии 
пластичности приводит к тому, что не только перепад давления между внешней и внутренней границами 
сферической оболочки, но и значения давлений на этих границах могут изменяться в ограниченном диапазоне. 
В данной постановке задачи, когда имеет место только тепловое воздействие, полый шар полностью не 
переходит в пластическое состояние. Результаты исследования позволяют прогнозировать поведение объекта 
(полого шара), испытывающего центрально-симметричные распределенные силовые и тепловые внешние 
воздействия. 


Ключевые слова: полый шар, толстостенная сферическая оболочка, термоупругопластическое состояние, 
эквивалентное напряжение, ассоциированный закон пластического деформирования, годограф напряжений, 
параметры управления поведением модели. 
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Введение. Решение задачи о толстостенной сферической оболочке, испытывающей разные внешние 
воздействия, приводится в монографиях [1,2] и ряде научных статей по теории упругости, пластичности, 
термоупругопластичности [3—9]. Обычно рассматривается случай, когда процесс нагружения является простым. 

Задача о толстостенной сферической оболочке является одной из простейших упругопластических задач, 
когда поля внешних воздействий и внутренних параметров обладают центральной симметрией. В силу 
центральной симметрии в пластической области выполняется режим полной пластичности. Для идеального 
пластического тела задача статически определимая, что позволяет ее решить для любых условий пластичности. 
В монографии [1] приведено наиболее полное решение задачи о сферической оболочке, когда не учитывается 
зависимость предела пластичности от температуры, а условие пластичности не зависит от первого инварианта 
тензора напряжений и знака третьего инварианта девиатора напряжений. Рассматриваются случаи только 
теплового и комбинированного нагружения, когда задана температура на стенках сферы, давление на 
внутренние стенки, а на внешней стенке давление отсутствует. В работах [10—13] выполнено исследование 
термоупругопластического состояния разных объектов, а в [3-9] рассмотрен процесс теплового нагружения и 
разгрузки свободного от внешних усилий шара и полого шара для условия Треска с учетом зависимости 
предела пластичности от температуры. Решение этой и аналогичных задач представляет интерес, поскольку 
можно получить аналитическое или частично аналитическое решение для разных математических моделей. 
Аналитическое решение можно получить при выборе кусочно-линейных функций пластичности [11, 13]. 
Математическое моделирование объектов позволяет прогнозировать их состояние и поведение в зависимости 
от значений исходных параметров [ 14, 15]. 


Механика 


23 


ЬИр://уез-допза.га 


А4уапсеа Епотеетт? Кезеагсй 2021. Т. 21, № 1. С. 22-31. 155М 2687—1653 





Материалы и методы. Постановка задачи. Рассматривается задача о толстостенной сферической 
оболочке (полом шаре), испытывающей центрально симметричные внешние воздействия: давление р, на 


внешнюю стенку при р=Ь и давление р, на внутреннюю стенку при р=а. Также рассматривается тепловое 
воздействие на шар: на границе р =а поддерживается температура Г, , на границе р=Б — температура Т,. 
Предполагается, что шар проявляет упругие и пластические свойства. Искомыми параметрами состояния в 
каждой точке шара являются компоненты тензора напряжений, компоненты тензоров деформаций и векторы 
перемещений. В области упругого состояния упругие деформации являются полными (остаточные деформации 
отсутствуют). 

Основные соотношения. Все соотношения приведены к безразмерному виду. В качестве масштаба 
длины выбирается внешний радиус шара Б. Все величины, имеющие размерность напряжений, отнесены к 
пределу пластичности при одноосном растяжении К. Масштабная единица для температуры — 1 С. 

В силу указанной симметрии внешних воздействий, в сферической системе координат р,0,ф матрицы 


компонент тензора напряжений и деформаций будут иметь вид: 


о 0 0 = 0 0 
(0)= 0 « 0, (=)= 0 в 0 
О 0 о ОО = 


ф ф 


В данном случае выполняются равенства 0, =0,, & =& 


ф? ф* 


Если функции пластичности не зависят от первого инварианта тензора напряжений и знака третьего 
инварианта девиатора напряжений, то при решении задачи о шаре функции пластичности будут приводиться к 
виду: 

Л=°,-©, ЕК. (1) 


Рассмотрим условие пластичности: 


(а, — 67)" +2(0, —85:)")" | 


о ,о,, 57,5 )= 
Л( р 0 р —. с+1(1+ о)" 


(2) 
. 1((, -с,-5(= — =’ |" +а(с, -в,-5(=’ —=)")" 
с+п + а)" 


— компоненты тензора пластических деформаций; ТГ — температура. 


= АСТ), 


Р 
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где 5, , & 
Когда параметры имеют значения: с =0, ц=1, 6 =0, а = 0, т=1 и=0,А=А, из условия (2) следует 
условие (1). На рис. | в плоскости о,,0, представлены кривые пластичности, определяемые по формуле (2) для 


разных значений числовых коэффициентов в функции пластичности. 


<] 
о, =—,ед 
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Рис. 1. Кривые пластичности: а — для параметров: с=0,2; т=2;ч= 0,5; д=0; т=3; А =1; сплошная линия 


© = 0,5; пунктирная линия © = 0; 6 — для параметров: = 0; с =0;0=0; А =1 


Артемов М. А. и др. Задача о толстостенной сферической оболочке 





Результаты, представленные на рис. |, показывают, что при учете первого инварианта тензора 
напряжений радиальное и окружное напряжения, когда точка шара находится в упругом состоянии, могут 
изменяться в ограниченном диапазоне. Соответственно давление на границы шара также должно быть 
ограниченным. Когда первый инвариант в условии пластичности не учитывается, упругое состояние возможно 
для любого значения давления на границах шара, но ограничен перепад давления Др=р, -р,. Учет знака 
третьего инварианта девиатора напряжений, как отмечалось выше, влияет на значения пределов пластичности. 

Если значения параметров состояния ©,,0, определяют точку области, ограниченную кривой 


пластичности, принимается, что определяющими уравнениями, связывающими напряжения и деформации, 
являются соотношения закона Дюамеля-Неймана [1, 2]: 


Ев, = (1-у)0,-уб,+ЕаТ, Её, =в, -2ув, + ЕОТ, (3) 
где модуль Юнга Ё и коэффициент Пуассона У — константы. 
Если параметры состояния о,,с, определяют точки на кривой пластичности, то принимается аддитивное 
представление полных деформаций через обратимые и необратимые деформации: 
де Р зе р 
5 = +5, &=5 +5, . (4) 


Полные деформации определяются через перемещения по формулам: 





и Чи 
# ==, & =—. (5) 
`а 
р р 
Полные деформации связаны условием совместности деформаций: 
ЧЕ 
го +&-& =0. (6) 
р 
Чи 
Приращения необратимых деформаций связаны с напряжениями нормальным законом: 
ЧЕ? Ч? 
ое (7) 


д//до, 0//до,` 


Соотношение (7) при выборе нелинейных функций пластичности в общем случае неинтегрируемое [16]. 
В квазистатическом приближении напряжения должны удовлетворять уравнению равновесия: 





Чо’, 
р +2(6,-6,) =0. (8) 
ар 
Эквивалентное напряжение. Эквивалентное напряжение — это выпуклые изотропные скалярные 


функции тензора напряжений. В частных случаях термин «эквивалентное напряжение» является синонимом 
других терминов, например, «интенсивность напряжений» [17]. В настоящей работе эквивалентное напряжение 
совпадает с функцией пластичности. В этом случае на упругопластической границе эквивалентное напряжение 
не будет иметь разрыва. 

Поле температур. Поле температур в шаре находится из решения краевой задачи [ | |: 


а О 
ар ар (9) 


и, 
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Решение задачи (9) представим в виде: 





г-т+-4АР В 1 АГАТ т. (10) 
(Б-а)\р 
Упругая область. В области упругого состояния полого шара формулы для напряжений имеют вид: 
В Ест 
с. а, в Е Е 
р р 2р 2р И-у)(6-а) 


Пластическая область. Выберем условия (1). Рассмотрим случай только теплового воздействия (10). 


Тогда пластическая область будет зарождаться на внутренней границе оболочки при условии [1]: $ 
в а дв _ 249 чай г 

(1—у) Ба+2Ь) з 

Обозначим через с — радиус упругопластической границы р=с. В процессе нагружения, когда 2 


В >В, пластическая область а <р<с, увеличивается. При выборе условия (1) напряжения в области а<р<с 


вычисляются по формулам: 25 
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о, =2кАШш(р/ а), ©,’ =’ +кА, к, =511(0,-0.), 
где к, =51(0, -в,)|.. .Если АТ >0,то к, =-1, если АТ <0,то к, =+1. 
Если область с <р<Ь остается упругой, то величины А, В и радиус упругопластической границы 
с, определяются из условий непрерывности напряжений на упругопластической границе и граничного условия 
о |= О. Так, если А и В определяются только из условий непрерывности напряжений на 


упругопластической границе, то имеют место выражения: 


. 
А=2к к ь [2+1 а ве ки: (11) 
Нл эл 56 33 





Уравнение для вычисления с, будет иметь вид: 


3 2 
2к^ ь (2+5 + а вв Л=0. (12) 
а) 3 3’ 3< 3’ Ь 


Если А и В определяются из условий непрерывности напряжений на упругопластической границе и условия 








с |, =0, то имеют место выражения: 


А= 











Е т м и, (13) 
3] 36 33 


Выбор формул (11) или (13) влияет на шаги алгоритма решения задачи, но не влияет на окончательные 
результаты. 
На границе р =ь будет зарождаться вторая пластическая область, если выполняется условие: 
(9,-6,) | ,=к,К, к, =—к,. (14) 
Для определения значения ДТ =ДТ, когда выполняется условие (14), необходимо решить совместно систему 
уравнений (12), (14). Поскольку параметр В в уравнения (13) и (1) входит линейно, то можно получить 
отдельно уравнение для определения радиуса упругопластической границы: 


дк А ОХе 9 —0 
а З(Б+с,)Бс, 


2 


(15) 


а также формулу для вычисления параметра В: 
в о 
° (5? —с”)аБ 


В процессе дальнейшего нагружения, когда выполняется неравенство В > В,, область шара с, <р<Ь 


(16) 


переходит в пластическое состояние на границе р=Ь. 
Полый шар при тепловом воздействии. Учет упрочнения материала. Рассмотрим случай, когда 
условие пластичности имеет вид: 
АО О Е, в Ле, = 56, =. ЛЕК, (17) 
Если до нагружения в шаре остаточные деформации отсутствуют, то в результате теплового нагрева 
пластическая зона будет зарождаться на внутренней границе р =а ‚ когда выполняются условия (15), (16). При 


дальнейшем нагружении образуется пластическая область а<р<с. Для нахождения напряжений в этой 


области надо из системы уравнений (3—6), (8), (17) получить соответствующие уравнения: 








‚ 4`в, 4, 6 баБЕоАТ 
р + 4р |= 0, 
Ар ар 1+30(1-уУ) р(Б-а) (18) 
рас, 
о, ==—+0.. 
2 ар 
Решение системы (18) запишем в виде: 
и ыы Уре сРвивАт |5 т 
+30(1-у —а 
р р (19) 





ИЕ И +2 Ш А + ы ВС 
1+361-у) 2р(Ь-а) 2р’ 
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Величины С,, С, ‚, входящие в формулы (19), определяются из граничного условия ©, | ,=0 и условий 
отсутствия пластических деформаций на упругопластической границе р=с: 


С =2к (К-К,)с' -ЗЕ а М, с», 


Я) — ° Е . 
С = кА па №), В№ 3“ р (20) 
: | За` а 


— барАТ 
° (+35а-у)6-а). 
В результате подстановки (20) в (19) получим: 
ЗбарЕалТ | 11 | _ 


1 
о ——_—_—_—_—__—_—_ НН ВНЕ 
` 1+360@-\) р-а а р 


2к Ас’ бабс,ЕаАТ 1 1 кет 1 
_ и |, 
3 р-а а р’ 3 Цар 


Ка ] | 


[кор 
а 





с, Е ВЕ кА+2кАшШР+ 
1+360(1-\) а р-а 


к Ас’ бабе ЕаАТ || 2 1 кА [2 1 
_ | ее 
3 2(Ь-а) а р’ за р 


Из решения упругой задачи следует: если ДТ > 0,то к =-1[; если ДТ <0,то к =+1. 


а 2р 





Учег пластической сжимаемости. Рассмотрим случай, когда функция пластичности является 
линейной относительно компонентов тензора напряжений: 
6(26. +6,)+к(5,-с,)+1(0, -о,)=&(-ХГ), 1 
к =51и(6, —С,). 
Условие пластичности (21) можно представить в виде: 
ас, +Во, =А(1-ХГ), 
&=2с+к-—\1, 
В=с-к-+1. 
С учетом введенных обозначений, для определения напряжений в пластической области получаем 
задачу: 
ос, +В, = (1-Х), 


Чо, 
Ро 620, (22) 


и 7 


Решение задачи (22) имеет вид: 


| (М -П 2МЕ |< ] г 1 2м | 
о =| р | № " + - К, 
Г " а-+В а(&-+2В) Др о+В р(а+28) 


В (М. -П  2М® а)“ [М- М, р 
= р | — + А, 
" (9) — о-в а(&-+2В) Др а+В  р(а+28))‘ 

адТ м аБАТ 
Ь-а’ р-а 


Для получения правильного результата из (20), когда, например, @+В=0, надо выполнить 





(23) 











где введены обозначения: М = Т, — ‚› М, =ХМ, М =ХМ. 


предельный переход при решении (22). Более просто получить правильный результат можно непосредственно в 
(23), при этом учесть, чтобы выполнялось условие @+[р}=0. В этом случае с = 0, поэтому получим: 


2 
с, 52% д-р [1-1 |2 | 
(9 а (ра 


с, А а-моанор+х, [1-2 =. 
(9) а р а 


В процессе нагружения пластическая зона зарождается на границе р =а ‚ когда выполняется условие: 
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_ @в+к-п(а- 25)БЕ а 
21-—у/(а? +аБ +?) } 
Результаты исследования. На рис.2 приведены графики напряжений и годографов вектора 


напряжений, когда область шара, соответствующая условию а <р < с ‚ находится в пластическом состоянии, а 


Ага -я в 


области шара, соответствующая условию с <р<ьЬ, — в упругом состоянии. 





ВЕТ 5 
а) 6) 
Рис. 2. Графики напряжений (а) и годографов вектора напряжений (6) 


при значениях параметров: А =; у = 0,3; а = 0,5; Б =ВАТ =170; с, =0,57 


На рис. 3 приведены графики напряжений и годографов вектора напряжений, когда области шара, 


соответствующие условиям а<р<с и с, <р<Ь, находятся в пластическом состоянии, а область шара, 


соответствующая условию с <р< с, , — в упругом состоянии. 
о {е] 
га , ед 9 = ра ь ед 





а) 0) 


Рис. 3. Графики напряжений (а) и годографов вектора напряжений (6) 


при значениях параметров: У = 0,3; а = 0,5; ДГ = 270; с, = 0,62; с, = 0,88 


На рис. 4 приведены графики напряжений и годографов вектора напряжений, когда область шара, 


соответствующая условию а<р<с ‚, находится в пластическом состоянии, а область шара, соответствующая 


условию с <р<Ь, — в упругом состоянии. 
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а) 0) 
Рис. 4. Графики напряжений (а) и годографов вектора напряжений (6) 
при значениях параметров: У = 0,3; а = 0,5; АТ = 215; с =0,58 
На рис. 5 приведены графики напряжений и годографов вектора напряжений, когда области шара, 
соответствующие условиям а<р<с и с, <р<Ь, находятся в пластическом состоянии, а область шара, 
соответствующая условию с <р<с, , — в упругом состоянии. 


7. ед о =—,ед 





а) 0) 


Рис. 5. Графики напряжений (а) и годографов вектора напряжений (6) при значениях параметров: 


у= 0,3; а = 0,5; АГ = 270; с, =0,61; с, =0,86 


На рис. 6 приведены графики напряжений и годографов вектора напряжений, когда области шара, 


соответствующие условиям а<р<с и с, <р<Ь, находятся в пластическом состоянии, а область шара, 


соответствующая условию с <р<с, — в упругом состоянии. 
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а) 0) 


Рис. 6. Графики напряжений (а) и годографов вектора напряжений (6) 


при значениях параметров: У = 0,3; а = 0,5; АТ = 79; с = 0,1; ч = 0,1; Еа = 0,012; ух = 0,0017 


Обсуждение и заключения. Результаты вычислений показывают, что в данной постановке задачи, 
когда имеет место только тепловое воздействие, полый шар полностью не переходит в пластическое состояние 
(рис. 2—6). Упрочнение приводит к увеличению эквивалентного напряжения в пластической области и 
уменьшению радиуса упругопластической границы (рис. 4, 5). Упругая область в процессе нагружения не 
может полностью исчезнуть. Пластическая сжимаемость и зависимость предела пластичности от температуры 
оказывают существенное влияние на напряженное состояние полого шара (рис. 6). 
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